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A new numerical method, recently developed to study ground states of 
the Falicov-Kimball model (FKM) , is used to examine the effects of correlated 
Q ■ hopping on the ground-state properties of this model in two dimensions. It is 

O , shown that the ground-state phase diagram as well as the picture of metal- 

insulator transitions found for the conventional FKM (without correlated hop- 
^N ■ ping) are strongly changed when the correlated hopping term is added. The 

ly-v . effect of correlated hopping is so strong that it can induce the insulator-metal 

^^ I transition, even in the strong-coupling limit, where the ground states of the 

J^ ' conventional FKM are insulating for all /-electron densities. 
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1 Introduction 

Since its introduction in 1969, the FKM |]I[ has become an important standard 
model for a description of correlated fermions on a lattice. The model was originally 
proposed to describe metal-insulator transitions and has since been investigated in 
connection with a variety of problems such as binary alloys p|, the formation of ionic 
crystals 0, and ordering in mixed- valence systems |§[. It is the latter language we 
shall use here, considering a system of localized / electrons and itinerant d electrons 
coupled via the on-site Coulomb interaction U . The Hamiltonian of the spinless 
FKM is 

H = Y, U.dtdj + f/ ^ f^f^dtd, + EfY. fth. (1) 

ij i i 

were f^, fi are the creation and annihilation operators for an electron in the local- 
ized state at lattice site i with binding energy Ef and df, di are the creation and 
annihilation operators for an electron in the conduction band. The conduction band 
is generated by the hopping matrix elements tij, which describe intersite transitions 
between the sites i and j. Usually it is assumed that t^ = —t ii i and j are nearest 
neighbors and t^j = otherwise (the conventional FKM), however, in what follows 
we consider a much more realistic type of hopping, so for the moment we leave it to 
be arbitrary. 

Recent theoretical works based on exact numerical and analytical calculations 
showed that the FKM, in spite of its relative simplicity, can yield the correct physics 
for describing of such fundamental phenomena as valence- change transitions, metal- 
insulator transitions, crystallization, charge ordering, etc. For example, it was found 
that the spinless FKM, in the pressure induced case, can describe both types of 
intermediate-valence transitions observed experimentally in rare-earth compounds: 
a discontinuous insulator-insulator transition for sufficiently strong interactions [Q 
and a discontinuous insulator-metal transition for weak interactions . In addition, 
at nonzero temperatures this model is able to provide the qualitative explanation 
for anomalous large values of the specific heat coefficient and for extremely large 
changes of electrical conductivity [^ found in some intermediate-valence compounds 
(e.g., in SmBg). Moreover, very recently the spin-one-half version of the FKM has 



been used to describe a discontinuous intermediate-valence transition (accompanied 
by a discontinuous insulator-metal transition) in SmS ^ as well as for a description of 
an anomalous magnetic response of the Yb-based valence-fluctuating compounds 0. 
On the other hand it should be noted that the model Hamiltonian (1) neglects 
all nonlocal interaction terms, and thus it is questionable whether above mentioned 
results persist also in more realistic situations when nonlocal interactions will be 
turned on. An important nonlocal interaction term obviously absent in the con- 
ventional FKM is the term of correlated hopping, in which the d-electron hopping 
amplitudes tij between neighboring lattice sites i and j depend explicitly on the 
occupancy {ftfi) of the /-electron orbitals. To examine effects of this term on 
ground-state properties of the two-dimensional FKM we choose the following form 
for the nearest-neighbor matrix elements 

Uj=U,+^^,{ftU + ftfi), (2) 

which represent a much more realistic type of electron hopping than the conventional 
hopping. 

Thus the spinless FKM in which the effects of correlated hopping are included 
can be written as 

h=y: u.did, + E t^(/+/. + f;f))dtd, + uY.ftf^dtd. + EfY. ftu (3) 

The first term of (3) is the kinetic energy corresponding to quantum mechanical 
hopping of the itinerant d-electrons between the nearest-neighbor sites i and j. The 
second term is just the correlated hopping term discussed above. The third term 
is the on-site Coulomb interaction between the d-hand electrons with density n^ = 
X Si dfdi and the localized /-electrons with density rif = j^Y^i ft fi^ where L is the 
number of lattice sites. The last term stands for the localized / electrons whose 
sharp energy level is Ef. 

Since in this spinless version of the FKM without hybridization the /-electron 
occupation number ft fi of each site i commutes with the Hamiltonian (3), the /- 
electron occupation number is a good quantum number, taking only two values: 
Wi = 1 oi 0, according to whether or not the site i is occupied by the localized / 
electron. Therefore the Hamiltonian (3) can be written as 



H= ^h,,{w)dtd, + Ef^w., (4) 



where hiAw) = tiAw) + UwiSa and 
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Uj (w) =tij + t[j (wi + Wj). ( 5 ) 

Thus for a given /-electron configuration w = {wi,W2 ■ ■ -Wl} the Hamiltonian 
(4) is the second-quantized version of the single-particle Hamiltonian h{w), so the 
investigation of the model (4) is reduced to the investigation of the spectrum of h for 
different configurations of / electrons. Since the d electrons do not interact among 
themselves, the numerical calculations should precede directly in the following steps, 
(i) Having w = {wi,W2 ■ ■ -Wl}, U, Ef and the nearest-neighbor hopping amplitudes 
t and t' fixed, (in the following t = —1 and all energies are measured in units of t) find 
all eigenvalues A^ of h{w). (ii) For a given Nf = J^i^a determine the ground-state 
energy E{w, U, Ef) = J2k=i ^ ^k + EfNj of a particular /-electron configuration w by 
filling in the lowest N^ = N — Nf one-electron levels (here we consider only the case 
Nf + Na = L, which is the point of the special interest for valence and metal-insulator 
transitions caused by promotion of electrons from localized / orbitals (/" — > /"~^) to 
the conduction band states), (iii) Find the w^ for which E{w, U, Ef) has a minimum. 
Repeating this procedure for different values of U,t' and Ef, one can study directly 
the ground-state phase diagram and valence transitions (a dependence of the /- 
electron occupation number on the /-level position Ef) in the FKM with correlated 
hopping. 

A direct application of this method has been used successfully in our previous 
papers p, [^ for a description of ground-state properties of the one-dimensional FKM 
model without correlated hopping {t' = 0). It was shown that finite-size effects 
are negligible for a wide range of the model parameters (e.g., strong interactions) 
and thus results obtained on relatively small clusters {L < 30) can be satisfactory 
extrapolated to the thermodynamics limit {L -^ oo). Using this method we have 
described satisfactory the strong-coupling phase diagram as well as the picture of 
valence and metal-insulator transitions in the one-dimensional spinless FKM [^ with 
t' = 0. It was found that for sufficiently large U the spinless FKM undergoes only a 



few discrete intermediate-valence transitions. These intermediate-valence transitions 
are insulator-insulator transitions, since they are realized between the insulating 
ground states corresponding to the most homogeneous configurations, which are the 



ground states in this region |rU|. Thus, there are no insulator- metal transitions in the 
1-d conventional FKM for strong interactions. In the next paper |]TT| we have shown 
that this picture of valence and metal-insulator transitions is dramatically changed 
if the term of correlated hopping is included. One of the most important results 
found for the one-dimensional FKM with correlated hopping was that the correlated 
hopping can induce the insulator-metal transition, even in the half-filled band case 
nd + nf = 1. In this paper we try to show that the same result holds also for the two 
dimensional case. Similar calculations are performed also away from the half-filled 
band case with the goal to examine possibilities for metal-insulator transitions in 
the strong-coupling limit. Another inspiration for performing these calculations was 
the recent paper of Wojtkiewicz and Lemanski [|I2|, where the authors studied two- 
dimensional FKM with correlated hopping using the combination of the perturbation 
expansion (up to the second order) and the method of restricted phase diagrams. 
They found that only a few phases form the ground-state phase diagram of the 
model in the strong coupling limit. For example, the ground state of the model for 
Ef = is the chessboard charge- density- wave (CDW) phase for all < t' < 1. Here 
we show that some other configurations, (e.g., the segregated configuration) can also 
be the ground states of the FKM at Ef = 0, thereby the ground-state phase diagram 
as well as the picture of metal-insulator transitions are strongly changed. 

2 The method 

Since the number of configurations that should be examined to obtain the ground 
state energy of the FKM grows exponentially with the system size, a direct applica- 
tion of the exact-diagonalization method described above is restricted to clusters up 
to 30 sites. In our previous papers we showed that clusters of this size are sufficient 
to suppress finite-size effects in one-dimension [^ ^, however, to obtain trustworthy 
results on the ground-state energy of the model in two dimensions one has to examine 
much larger clusters {L ~ 100). Unfortunately, the clusters with L > 30 are beyond 



the reach of present day computers within exact diagonahzations, and thus the only 
way is to compute the ground-state properties of the model by an approximate but 
well controlled method. Here we use the simple method based on a modification of 
the exact-diagonalization procedure described above. The method consists of follow- 
ing steps: (i) Chose a trial configuration w = {wi,W2 . . .wl}. (ii) Having w, U and 
Ef fixed, find all eigenvalues A^ of h{w) = T + UW. (iii) For a given Nf = J2i Wi de- 
termine the ground-state energy E{w) = J2k=i * ^k + EfNf of a particular /-electron 
configuration w by filling in the lowest Nd = L ~ Nf one-electron levels, (iv) Gener- 
ate a new configuration w' by moving a randomly chosen electron to a new position 
which is chosen also at random, (v) Calculate the ground-state energy E{w'). If 
E{w') < E{w) the new configuration is accepted, otherwise w' is rejected. Then the 
steps (ii)-(v) are repeated until the convergence (for given parameters of the model) 
is reached. Of course, one can move instead of one electron (in step (iv)) two or 
more / electrons, thereby the convergence of method can be improved. Indeed, tests 
that we have performed for a wide range of the model parameters showed that the 
latter implementation of method, in which A^o > 1 electrons (A'"o should be chosen 
at random) are moved to new positions, overcomes better the local minima of the 
ground-state energy. This also improves the accuracy of method. 



This method was first used in our recent paper [|13l to study the ground-state 
properties of the one and two-dimensional FKM without correlated hopping. It 
was found that on small and intermediate clusters {L ~ 30) the method is able to 
reproduce exactly the ground states of the conventional FKM, even after relative 
small number of iterations (typically 10000 per site). For such clusters the method 
is only rarely stopped at the local minimum. Of course, with increasing L the 
problem of local minima appears often. Fortunately, it can be considerably reduced 
by more efficient algorithm (one is discussed above) or by increasing the number 
of iterations. The latter case imposes, however, severe restrictions on the size of 
clusters than can be studied with this method {L ~ 100, for 10^ iterations per 
site) . To verify the convergence of this method for the two-dimensional FKM with 
correlated hopping we have performed the same calculations on the cluster of 4 x 4 
sites, where ground states can be obtained also within the exact diagonalization 
calculations. Numerical results obtained for a wide range of the model parameters 
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{t' = —1, —0.8, . . .1,U = 0, 0.1 . . . 10) shoved that the exact ground-states can be 
again reproduced after ~ 10000 iterations per site. 

3 Results and discussion 

The most interesting question that arises for the FKM with correlated hopping is 
whether the correlated hopping term can change the ground-state phase diagram 
and the picture of valence and metal-insulator transitions found for the conventional 
FKM (t = — 1 and t' = 0). The nature of the ground state, its energetic and struc- 
tural properties, and the correlation-induced metal-insulator transitions are subjects 
of special interest. For the conventional FKM these problems are well understood 
at least in the symmetric case {Ef = 0,nf = n^ = 1/2). In this case the localized 
/-electrons fill up one of two sublattices of the hypercubic lattice (the charge-density- 
wave state) and the corresponding ground state is insulating for all f/ > 0. Thus, for 
the finite interaction strength there is no correlation-induced metal-insulator transi- 
tion in the symmetric case. 

One can expect, on the base of simple arguments, that the ground-state phase 
diagram of the FKM with correlated hopping will be fully different from one dis- 
cussed above for the conventional FKM. Indeed, the following selection of hopping 
matrix amplitudes t = — 1 and t' > may favor the segregated configuration since 
the itinerant d electrons have the lower kinetic energy in this state. This mechanism 
could lead, for example, to the instability of the CDW state that is the ground state 
for t' = 0, and thereby to a metal-insulator transition, even in the symmetric case. 
To examine possibilities for such a transition in two dimensions we have performed 
an exhaustive study of the model on 6 x 6 and 8x8 clusters (with periodic bound- 
ary conditions) for a wide range of parameters t' and U. The results of numerical 
calculations are summarized in Fig. 1 in the form of the t'-U phase diagram. In 
addition to the CDW state wi that is the ground state at t' = for all nonzero U we 
found two new phases that can be the ground states of the model, and namely, the 
configuration with alternating lines of occupied and unoccupied sites W2 and the seg- 
regated configuration ^3 (see Fig. 2). Thus at nonzero t' the CDW state wi becomes 
unstable against the transition to W2 and w^. The transition from wi to W2, as well 



as from Wi to w^ is the insulator-metal transition since the configuration wi has the 
finite gap (~ U) at the Fermi energy ^^ for all nonzero values of U, while both W2 
and W3 are metallic in the corresponding regions of stability. Thus we arrive to the 
very important conclusion, and namely, that the correlated hopping term can induce 
the insulator-metal transition, even in the half-filled band case. Another important 
result, confirming the crucial role of the correlated hopping term can be seen from 
Fig. 1, where the comprehensive phase diagram of the two-dimensional FKM with 
correlated hopping is presented. It is seen that the correlated hopping can destroy 
the CDW state, even at large values of the Coulomb interaction ([/ ~ 7). This is 
an unexpected result since recent results of Wojtkiewicz and Lemanski |jl2| based 
on the combination of the perturbation expansion (up to the second order) and the 
method of restricted phase diagrams predicted that the ground state of the model at 
Ef = and U large is the CDW state for all values of < t' < 1. This discrepancy 
is probably due to the fact that the authors examined (as possible ground states) 
only a restricted set of configurations (consisting of all periodic configurations having 
elementary cells up to 12 sites), and the segregated configuration (that should be the 
ground state in this region) does not belong to this set. Another possible explanation 
of this discrepancy is that the second-order perturbation expansion used by authors 
is insufficient to describe correctly the ground-state properties of the model in this 
region (f/ ~ 7). 

The fact that the correlated hopping can induce metal-insulator transitions indi- 
cates that the picture of valence and metal-insulator transitions found in our previous 



papers within the conventional FKM [g, 0, 13| should be dramatically changed if finite 



values of t' will be considered. The largest changes are expected in the strong cou- 
pling limit (f/ > 4), where all ground states of the conventional FKM are insulating 
for both 1-d and 2-d case 0, 10, 13], while the numerical results obtained for nonzero 
t' show on the existence of the metallic phase, at least for nj = 1/2. We suppose that 
this important result is not restricted to the half-filled band case only, but persists 
also for /-electron densities away from this point. To verify this conjecture we have 
performed an exhaustive study of the model for nj = 1/4 on 6 x 6, 8 x 8 and 12 x 12 
clusters. Our numerical calculations showed that the phase diagram of the model at 
Uf = 1/4 is separated into two distinct regions. In the first region {U < 2) the phase 



diagram has a complex structure with the ground state apparently changing point 
by point at every value of the correlated hopping amplitude t' for fixed interaction 
strength. Unfortunately, the structure of the phase diagram in this region strongly 
depends on the size of cluster and thus we were not able to extrapolate satisfactory 
these results to the thermodynamic limit L -^ oo. Contrary to this case, the phase 
diagram exhibits a very simple structure (see inset in Fig. 1) in the opposite limit 
{U > 2). In this region only two configurations are the ground states of the FKM 
with correlated hopping, and namely the segregated configuration and the configu- 
ration W4 (see Fig. 2) that has been proven to be ground state of the conventional 



FKM for large U (see Ref. [|13], |T5|, |T6[). Since the configuration Ws is metallic and w^^ 
insulating we have the correlated hopping induced metal-insulator transitions also 
at Hf = 1/4. The metallic phase is stable up to f/ ~ 7 and this again confirms our 
conjecture that the comprehensive picture of metal-insulator transitions in the FKM 
with correlated hopping will be fully different from one found for the conventional 
FKM, especially for U large. To complete this picture one has to perform similar 
calculations for all /-electron densities what is a cumbersome computational task, 
even on 8 x 8 cluster. The work on this subject is currently in progress. 

In summary, the effects of correlated hopping on the ground-state properties of 
the FKM in two dimensions have been studied. It was shown that the ground- 
state phase diagram as well as the picture of metal-insulator transitions found for 
the conventional FKM are strongly changed when the correlated hopping term is 
added. The effect of correlated hopping is so strong that it can induce the insulator- 
metal transition, even in the strong-coupling limit, where the ground states of the 
conventional FKM are insulating for all /-electron densities. 

This work was supported by the Slovak Grant Agency VEGA under grant No. 
2/7021/20. Numerical results were obtained using computational resources of the 
Computing Centre of the Slovak Academy of Sciences. 
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Figure Captions 

Fig. 1. t'-U phase diagram of the two-dimensional FKM with correlated hopping 
at half-filling [Ej = 0,nf = Ud = 0.5). Three different phases correspond to the 
CDW state wi, the configuration with alternating lines of occupied and unoccupied 
sites (^2), and the segregated configuration W3. The inset shows t' — U phase dia- 
gram for Hf = 1/4 and U > 2. Two different phases correspond to the segregated 
configuration and the configuration W4 that has been proven to be ground state of 
the conventional FKM for large U. 

Fig. 2. The ground-state configurations of the two dimensional FKM with corre- 
lated hopping for Uf = 1/2 (^1,^2 and W3) and Uf = 1/4 (^4). 
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